function [mean_X,std_X] = Calc_monthly_vol(stationary_prob,X)
% computes std(X(t)) along with intermediate statistics

    mean_X = sum(stationary_prob(:).*X(:));

    std_X = sqrt(sum(stationary_prob(:).*(X(:).^2)) - mean_X^2);

end